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Period-Magnitude relation of Mira-like variables in the Large 
Magellanic Cloud as a tool to understand circumstellar extinction 
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ABSTRACT 

Near- to mid-infrared period-magnitude relations and also the period-bolometric luminosity 
relation of OGLE-III Mira-like variables in the LMC are derived. The relations have a kink, 
and the period at which the break occurs is quantitatively obtained. There are many Mira-like 
variables whose fluxes at the optical and the near-infrared wavebands are fainter than the ones 
predicted by the period-magnitude relations. The deviation is due to the circumstellar extinc- 
tion, and the amount of the deviation is found to be strongly correlated with near-infrared 
colors. The empirical formulae relating the amount of the deviation and the near-infrared col- 
ors are derived. These relations are useful to accurately calculate the distances to the dusty 
Mira-like variables, because the dimmed fluxes due to the circumstellar extinction can be es- 
timated. In a manner analogous to the interstellar extinction law, the ratios of deviations at 
any two different wavebands are calculated. The ratios are found to change with the pulsation 
period, indicating that the dust properties are subject to change as Mira-like variables evolve. 
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1 INTRODUCTION 

One of the most astronomically profitable aspects to study radi- 
ally pulsating variable stars is that some groups of them have their 
individual period-luminosity relations that can be used as cosmic 
distance scales. Mira-like variables are one of such group of vari- 
able stars that evolved from low- to intermediate mass (~ 0.8 < 
M/M e < ~ 8) stars, and they ar e in the late stage of stellar evolu- 
tion (e.g. Jlben & Renzinil 19831) . Since the di scovery of the period- 
magn itude relation of Mira-like variables l lGlass & Llovd Evansl 
1 198 lh . many astronomers observed them to refine the relation and 
to study its metallicity dependency. From late 90's to early 00's, 
there was a major progress on the study of the period-luminosity 
relations of variable stars due to the advent of the optical large 
survey projects looking for the gravitational lensing events with 
robotic telescopes (e.g., MACHO, OGLE, EROS, MOA). The dis- 
covery of parallel sequences of period-magnitu de relations of red 
giant variables in the Large Magellanic Cloud jWood|[200fj|) is an 
especially important result, because it told us that each group of 
variables has different period-luminosity relations, and also, dif- 
ferent relation is assigned t o different pulsation mode within a 
group (see also, for example. iKiss & Bedding 2003, Ita et al. |2004 
Groenewegen et alj|2004l. iFraser et al J 1200a. iDerekas et alj|200q . 
Glass et alj2009LISoszvnski et alj|2009l, and lTabur et all201Ch . 



The aim of this paper is to study the near- to mid- 
infrared period-magnitude relations of Mira-like variables in the 
LMC, us ing the combined da ta of recently published OGLE- 
III data JSoszvnski etaP l2009h and the Spitzer SAGE catalog 
jMeixner et alj2006f) . Infrared data is imperative because Mira-like 
variables are usually associated with the circumstellar dust shell, 
which is formed by the mass-loss phenomenon. 



2 THE DATA 



Recently, the OGLE project released its third phase 

(OGLE-III; 12th lune, 2001 - 3rd May, 2009) survey data. They 
provided catalogs of several types of variable stars in the Mag- 
ellanic Clouds. Using their on-line databasj], stars with variable 
star type of "Mira" were chosen b y querying their cata log of 
long-period variables in the LMC Jsoszvnski etafl l2009h . The 
query yielded a catalo g of 1,663 Mira-like variables in the LMC. 
ISoszvrtski etai]j2009l) used the /-band pulsation amplitude to dis- 
tinguish Miras and Semi-regular variables. They defined Miras as 
having /-band pulsation amplitude greater than 0.8 mag. In this pa- 
per, we follow their classification and call these Mira-like variables 
as just Miras or OGLE-III Miras. The OGLE-III catalog lists not 
only coordinates, but also pulsation periods, amplitudes, surface 
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chemistries (Oxygen-rich / Carbon-rich) inferred by photometric 
colors, time averaged V- and /-band magnitudes, and other useful 
parameters. 



2.1 Cross-identification with existing catalogs 

The OGLE-III Mira catalog is cross-identified with the following 
existing catalogs using a positional tolerance of 3 arcsec. If more 
than one stars are present within the tolerance radius, the closest 
one is adopted and discard the others. The result is used for discus- 
sion in the rest of this paper. 

The Large Ma gellanic Cloud Photometric Survey catalog 



Table 1. The adopted zero magnitude flux densities and their corresponding 
reference wavelengths. 



,arge Ma ; 

dZaritskvet afll2004l) : The catalog lists U, B, V, and / stellar pho- 
tometry of the central 64 deg 2 area of the LMC. It must be noted 
that we preferentially used the time-averaged V— and /-band data 
in the OGLE-III survey catalog whenever available. 1,412 out of a 
total of 1,663 Miras have counterparts in this catalog. 

• The Two Micr on All Sky Survey (2MASS) catalog 
dSkrutskie etai]|2006l) : The catalog provides uniform /, H, and 
K s photometry for sources all over the sky. The 2MASS catalog 
is complete down to K s < 14.3 mag in the absence of confusiorj 2 ]. 
The 2MASS magnitudes are used for the bright sources without 
IRSF measurements (The IRSF survey did not detect bright sources 
due to the saturation limit of about 10 mag at A", band. See below). 
1,639 out of a total of 1,663 Miras have counterparts in this catalog. 

• The IRSF Magell anic Clouds Point Source Catalog (IRSF cat- 
alog; |Katoet^][2TO3): The IRSF catalog lists J,H, and K s pho- 
tometry of over 1.4xl0 7 sources in the central 40 deg 2 area of the 
LMC. Comp ared to the contempor ary DENIS dCioni et al]|2000h 
and 2MASS dSkrutskie etai1l2006l) catalogs, the IRSF catalog is 
more than two magnitudes deeper at K s band and about four times 
finer in spatial resolution. The IRSF system magnitudes were con- 
verted into the 2MASS system o nes b y using the conversion equa- 
tions given in IKato et al j d2007t) and IKucinskas et all d2008h . We 
preferentially used IRSF photometry whenever available. 1,557 out 
of a total of 1,663 Miras have counterparts in this catalog. 

• The Spitzer SAGE-LMC survey catalog dMeixner et al.l2006f) : 
The catalog lists near- ([3.6], [4.5], [5.8], and [8.0]) to mid-infrared 
([24] and [70]) photometry of sources in the central 49 deg 2 area 
of the LMC. Throughout this paper, the numbers bracketed by [ ] 
designates the data of the Spitzer catalog, for example, [3.6] indi- 
cates the photometry in the 3.6 pm band. The SAGE-LMC team re- 
cently made a final data delivery (DR30. They provided two types 
of the catalog. One is the "Catalog", which is a more highly reliable 
list. The other is the "Archive/Full", which is a more complete list. 
See the explanatory document "The SAGE Data Products Descrip- 
tion" prepared by Dr. M. Meixner for more detailed descriptions 
of the catalogs. We use the "Archive/Full" version for the cross- 
identification. 1,612 and 1,274 out of a total of 1,663 Miras have 
counterparts in the near- and mid-infrared SAGE-LMC catalog, re- 
spectively. 

• Optical carbon star catalog dKontizas et alj200lh : The catalog 
lists 7,760 carbon stars in the LMC that were identified by the opti- 
cal objective-prism spectroscopy survey. 163 out of a total of 1,663 
Miras have counterparts in this catalog. 

• Spitzer IRS-identifie d O-rich and C-ri c h evo lved stars 
dGroenewegen et aT]|2009h : iGroenewegen et al ] d2009l) uniformly 



http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec2_2.html 
visit http://ssc.spitzer.caltech.edu/legacy/sagehistory.html 
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reduced all Spitzer IRS data of evolved stars in the Magellanic 
Clouds taken so far. They calculated luminosities and mass-loss 
rates for 66 O-rich and 68 C-rich infrared stars in the LMC. 40 out 
of a total of 1,663 Miras have counterparts in their list. 

2.2 Corrections for the interstellar reddening 

In this paper, the optical (U, B, V, and /) and the near-infrared (J, H, 
and K s ) photometry are corrected for the interstella r redde ning 
based on the relations in ICardelli. Clayton. & Mafhisl dl989T) . as- 
suming A V /E(B - V) = 3.2. We uniformly adopt (Ay, A B , Ay, 
A,, Aj, A H , A Ks ) = (0.407, 0.345, 0.272, 0.159, 0.078, 0.048, 
0.032) mag, corresponding to the mea n foreground reddening of 
E(B- V) = 0.085 mag toward the LMC dLarsen. Clausen. & Storml 
I200CT) . Fluxes in any other longer wavelengths (i.e., photometry in 
the Spitzer wavebands) are not corrected for the interstellar extinc- 
tion, which we assume negligible. We also ignore the reddening 
inside the LMC. 

2.3 Calculation of the bolometric luminosity 

We calculate the bolometric luminosity of all OGLE-III Miras by 
using a cubic spline to interpolate the spectral energy distribution 
and integrate it from the shortest available wavelength to 8 or 24 
pm. We used zero-magnitude fluxes tabulated in Table[T]to convert 
magnitude into Jansky. Color-correction is not applied to the flux 
density, due to the lack of information of the incident spectrum. For 
some of the very red sources, the calculated luminosities can be 
underestimated to a large extent because the fluxes longward of 24 
pm are not included. Therefore the calculated luminosities should 
be only lower limits for the very red sources. Also, we have to note 
that the calculated luminosities are rather uncertain, because they 
are sensitive to changes of degree of interstellar reddening, color 
correction, and also to time variation of light, which is significant 
especially for Mira-like variables discussed here. 



2.4 Evaluation of the color-classified surface chemistry in the 
OGLE-III catalog 

ISoszvnski etall d2009l) used Wj vs. Wjk diagram to classify O-rich 
and C-rich surface chemistry of Mira variables in the LMC, where 
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Figure 1. The Wi vs. Wjk diagram of OGLE-III Miras in the LMC, where 
Wi and Wjk are Wesenheit indices (see text). The black and red small dots 
show the color-classified O-rich and C-rich chemistries (iSosz vnski et alj 
2009), respectively. Miras with spectroscopically known surface chemistry 
are shown in big pluses. Th e green and red pluses show optical and infrared 
C-rich stars identified by Kontizas et al. (2001) and Groenewegen et al. 
(2009), respectively. The blue plu ses indicate O-rich red supergiants identi- 
fied by Groenewegen et al. I200JJ). 



Wi and Wjk are Wesenheit indices defined as Wi = I — 1.55(V - /) 
and W JK = K s - 0. 686(7 - K s ), respectively. The validity of their 
classification should be checked before use, and its evaluation is 
done in the following way. W e make the same Wi vs. Wjk diagram 
as in lSoszvnski et all J2009I) . using only 1,663 OGLE-III Miras. 
The diagram is shown in figure [T] Then, those with spectrosc opi- 
cally known surface chem istry (based on iKontizas et alj|200ll and 
iGroenewegen et al]|2009h are highlighted to see their distribution 
on the employed Wesenheit index plane. It is clear that the distribu- 
tion of Miras with spectroscopically known chemistry is fairly re- 
produce d by the OGLE's color -classification criteria, that were de- 
fined bv lSoszvhski et all d2009h . We do not see any color-classified 
O-rich stars that correspond to the infrared O-rich red supergiants 
(blue pluses in Figure QJ. This is due to the lack of their V-band 
data, probably because of the circumstellar extinction. Presumably, 
they get fainter than the detection limit of the OGLE-III V-band 
observations. We use their color-classified chemistry throughout 
this paper unless otherwise noticed. 



3 RESULTS AND DISCUSSIONS 
3.1 Period-magnitude relations 

iGlass et all d2009l) used t he combined data of the MACHO survey 
(e.g.. lAlcock et alj|2000h and the pr evious version of the Spizer 
SAGE catalog dMeixner et alj|2006l) to discuss the mid-infrared 
period-magnitude relations (PMR) of variable red giants. They 
compared the PMRs of variable red giants in the LMC and the 
NGC6522, showing that there is little difference in their PMRs 
despite the differences in ages and metalli cities between the two 
galaxies. Very recently, iRiebel et al .1 d2010l) discussed the infrared 




Log 10 (Period) [days] 

Figure 3. The same as the Figure [2] but stars with spectroscopically 
known surface chemistry are highlighted. The pluses show op t ical ( green) 
and infrared (red) C-rich stars identified by Konti zas et alj bOOll) and 
Groenewegen et al 1 J2009T) . respectively. The blue cros ses indicates O-rich 
red supergiants identified bvlGroenewegen et al. 12009). 



PLRs and showed that the wavelength dependence of the slope of 
the period-luminosity relationship is different for different classes 
of variable red giants. 

Here we use the combined data of the OGLE-III Mira catalog 
and the new version of the Spitzer SAGE catalog (DR3), focus- 
ing on the PMRs of Miras in the LMC. In Figure [2] the PMRs of 
the OGLE-III Miras are shown. The abscissa is the common loga- 
rithm of the primary pulsation period in days. We use the primary 
period only. Note that the OGLE-III catalog provides not only the 
primary period, but also the secondary, and the tertiary periods. The 
left ordinate of each panel is in reddening-corrected apparent mag- 
nitude. Its corresponding waveband is indicated at the top left of 
each panel. The unit of the left ordinate can be scaled to the abso- 
lute one by subtracting a certain constant (i.e., distance modulus of 
the LMC). The colors of the marks in Figure [2] represent the dif- 
ference in the color-classified surface chemistry, as being black for 
O-rich, and red for C-rich Miras, respectively. Figure[3]is the same 
diagram as the Figure [2] but highlights Miras with spectroscopi- 
cally known surface chemistry. The plus es show optical (green ) and 
infrared (red) C-rich stars identified by I Kontizas et alj feOOll) and 
IGroenewegen et aT] d2009h . respecti vely. The blue crosses indic ates 
O-rich red supergiants identified bv lGroenewegen et alj J2009h . By 
comparing these two figures, one can understand where each of the 
three groups is located in the period-magnitude plane. 
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Log 10 (Period) [days] 

Figure 2. P eriod-magnitude relations for OGLE-HI Miras in the LMC. The color of the marks represents the differences in the color-classified surface 
chemistries I Soszvriski et al. 2009); black for O-rich and red for C-rich stars. The cyan vertical lines show the mean magnitudes of O-rich Miras in 0.05 mag 
bins and their corresponding standard deviations. The thick solid and thick dashed blue lines are least-square fit to the mean magnitudes (see text and table[3j. 
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Table 2. Results of Student's t-test to test a hypothesis that "the O-rich 
and C-rich Miras obey the same period-luminosity relation". The A and R 
means accept and reject the hypothesis, respectively. 
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There are several interesting features in Figures [2] and [3] It 
seems that Miras with shorter primary pulsation periods fall on 
a relatively tight PMR regardless of their surface chemistry. We 
conducted a two-tailed Student's t-test to test a hypothesis that 
"the O-rich and C-rich Miras obey the same period-luminosity 
relation". A 95 % confidence level was chosen. Specifically, the 
mean magnitudes and variances of O-rich and C-rich Miras in the 
log 10 P(days) = 0.05 bins are calculated and then compared. The 
results are summarized in Table|2] The test suggests that Miras with 
periods shorter than about 200 days (log 10 P ~ 2.3) obey the same 
PMR regardless of their surface chemistry and also of employed 
waveband. For the easy explanation, we call the PMRs of these 
shorter period Miras as standard PMRs, and use them for compari- 
son. 

The PMRs become complex for Miras with pulsation periods 
longer than about 320 days. C-rich Miras with longer periods locate 
below the extension of the standard PMRs in the optical and near- 
infrared wavebands, although they locate above the standard PMRs 
in the Spitzer wavebands. This can be understood as the extinc- 
tion by the circumstellar dust, and also as the intrinsic redness of 
C-rich Miras. On the other hand, O-rich Miras with longer periods 
locate above the reference PMR in all wavebands. Some O-rich Mi- 
ras with very long period (log 10 P > 3.0) locate below the reference 
PMR, due probably to the extinction by circumstellar dust, just like 
the C-rich Miras. Indeed, some o f them are spectroscopica lly con- 
firmed to have circumstellar dust dGroenewegen et al.l2009h . These 
dusty O-rich Miras with very long period may be counterparts to 
the Galactic OH/IR stars. 

To derive the PMRs, we concentrate on the O-rich Miras (i.e., 
black dots in Figure^. The scatter about the PMR is related to dif- 
ference in the color, error in the photometry, difference in the ob- 
served pulsation phase, and its intrinsic error. Mira-like variables 
discussed here are bright especially in the infrared, and their photo- 
metric errors should be small. Meanwhile, they are large-amplitude 
variables, and the difference in the observed pulsation phase (recall 
that we use single-epoch photometric data) can produce significant 
scatter. With the present data, we can not estimate the attributable 
fraction of scatter for each possible cause. Then we assume that 
the primary cause of the scatter about the PMR is due to the dif- 
ference in the observed pulsation phase, and derive the PMRs in 
the following way. As a first step, we calculate the mean apparent 
magnitude (< m, >) of O-rich Miras in the log 10 P (days) = 0.05 



bin and its corresponding standard deviation (o"„, ; ), where the suf- 
fix i stands for the ;-th bin. The cyan points with error bars in 
Figure [2] show the calculated < m, > and o~ m .. Because the cyan 
points obviously do not lie on a line, we fit two lines in the form 
of m A = a A log 10 P + b A to the cyan points by using least-square 
algorithm. We gave a relative statistical weight for each data point, 
which is proportional to 1/crjj,.. Apparently, the break seems to oc- 
cur at around log 10 P (days) = 2.6. In the following section, we 
estimate the break periods quantitatively. The calculated PMRs are 
shown as the blue thick solid and dashed lines in Figure[2] The val- 
ues of the a A and b A with their errors and the log 10 P range used 
for the fitting are summarized in Table [3] Also, we calculated the 
residuals r from the derived PMRs for all O-rich Miras, which is 
defined as r = a A x log P + b A - '^observed ■ ^ ne standard deviation 
of r is then calculated after 10 times iterations of 3 sigma clipping. 
The standard deviation of r and the number of O-rich Miras used 
to calculate it are also tabulated in Table [3] The PMRs of C-rich 
Miras are not derived because most of them (except optical ones 
with periods shorter than about 400 days) do not obey a relation in 
the optical and near-infrared wavebands. In addition, their PMRs in 
the Spitzer wavebands have very large scatter, probably due to the 
broad range in the amount of circumstellar extinction and/or in the 
amplitude of light variation. 



3.2 A kink in the period-magnitude relation 

It is clear that the cyan points do not fall on a single straight line. 
There is a break at a certain pulsation period, and the PMR becomes 
steep from there. This feature is seen regardless of the wavebands. 
We determined the break period by calculating a, = (< m i+i > 
- < nij >)/0.05, where the a f should remain constant (equals to the 
slope) for a straight line without a kink. Although the calculated a, 
is a bit noisy, we found a significant leap of a, at a certain pulsation 
period. The corresponding periods (log 10 P) for the leap are, 2.7 for 
<I>,J, H, K s , and m bol , 2.65 for [3.6], [4.5], [5.8], and [8.0], and 
2.6 for [24 ] , resp ectively. 

iFeastl < fl989h first pointed out the existence of O-rich Miras in 
the LMC with periods longer than 420 days (log 10 P ~ 2.6) that 
locate above the extrapolations of the period-A - magnitude rela- 
tion derived with Miras wit h periods shorter than 4 20 days. Then 
iHughes & Woodl dl990h and lWhitelock et al.1 d2003h increased the 
number of samples, showing that the period-ra bo i relation also have 
a break at around that period. Our results are very consistent with 
these previous works. A possible explanation for the excess lumi- 
nosity of O-rich Miras wit h longer periods is the "hot bottom b urn- 
ing" (HBB) process (e.g.. IWhitelock et alj|2003l lFeastll2009l) that 
occurs in stars with mass es greater than about 3 M e wit h the metal- 
licity of the LMC (e.g.. iGarcfa-Hernandez et al.|[2006h . Note that 
the lower mass li mit to activate the HBB depends on the m etallicity 
(see, for example. iBoothrovd. Sackmann. & Ahernll 1993b . Very in- 
terestingly, all but one spectroscopically-known infrared C-rich Mi- 
ras have periods longer than the break period. Pulsa tion period o f 
a Mira variable is related to its mass and radius (e.g., IWoodlll990h . 
Assuming that the radii of O-rich and C-rich Miras with about the 
same pulsation period are more or less the same, the longer period 
one should have larger mass. On the condition that the kink in the 
PMR is due to the HBB, one can speculate that the C-rich Miras 
with periods longer than the break period are undergoing the HBB 
process now. If so, the high mass-loss rate (> 10~ 5 M G /yr, see next 
section) of infrared C-rich Miras can be due to the excess luminos- 
ity emerged from the HBB process that should create additional ra- 
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Table 3. Period-luminosity relations for O-rich Mira variables in the LMC 
of the form mx = a a X log P + bn . 
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Figure 4. The correlation between the observed near-infrared color of all 
color-classified O-rich (black dots) and C-rich (red dots) Miras and the "de- 
viation" (see text). The solid line is the best-fit of a quadratic function in the 
form of d = ax 2 + bx + c, where the d and x are the deviation and the em- 
ployed near-infrared color in magnitude. The fitting results are summarized 
in Tableg] 



diation pressure to the dust grains. See, for example, lMarigdJ2007h 
for the theoretical view on the HBB process in C-rich AGB stars. 



Table 4. Color-deviation relations for C-rich Mira variables in the LMC of 
the form deviation [mag] = ax 2 + bx + c, where x is the employed color. 



deviation x a b c 

Hobs - Hpl H-K s 0.939±0.029 -0.038±0.077 -0.155±0.045 

H obs -H PL J-K s 0.205±0.010 -0.041±0.056 -0.301±0.071 

K Sob ,-K spL H-K s 0.924±0.029 -0.968±0.076 0.166±0.044 

K s °„j-K spL J-K s 0.213±0.009 -0.586±0.052 0.331±0.066 



3.3 Circumstellar extinction 

Figure [3] and Table [2] suggest that the optical O-rich and C-rich 
Miras fall on almost the same PMRs in optical and near-infrared 
wavebands, but infrared Miras do not. It is well known that the 
PMR is best defined at A"-band, where the atmospheric bands do 
not contribute much and the fluxes correspond to the photosphere. 
It seems that this explanation holds for the optical Miras only. 
The A^-band fluxes of the infrared Miras, however, do not rep- 
resent their photospheres any more be cause of the heavy circum- 
stellar extinction even in the K - band. I Wood | |2003|) studied MSX 
fegan. Van Dvk. & Pricd200lh sources in the LMC, and found that 
infrared Mira-like variables do not fall on the standard period-A - , 
magnitude r elation for M ira- like variables ( See, for example, se- 
quence C in lWoodl20r3(j| and llta et al1l2004l) , and they are fainter 
than the luminosity that the PMR predicts. Here we assume that 
the reason of the reduced optical and near-infrared flux of the in- 
frared Miras is due solely to the circumstellar extinction by their 
surrounding dust shell. We also assume that all of the infrared Mi- 
ras should have fallen on the extension of the optical and near- 
infrared PMR of optical Miras if they were without circumstellar 
dust. Under these assumptions, the deviation from the luminosity 
that the PMR predicts (i.e., m^ deviation = m^bserved " "Upmr) 
should be a measure of the strength of circumstella r extinction 
dMatsunaga. Fukushi. & Nakaddl2005l . Ilta et alj|2009h . It is inter- 
esting to note that there is no relevance between the deviation and 
pulsation period (see Figures [2] and [3} . It is also worth mentioning 
that the deviation is also independent of the pulsation amplitude at 
near-infrared bands. The interpretation of these observational facts 
will be made in the separate paper (Ita et al., in preparation). 

As stated above, the amount of deviation is independent of 
pulsation period and amplitude. However, it would be of great in- 
terest if we can correlate the deviation with the other observables 
(and especially distance-independent ones), because such relation 
would be useful, for example, to correctly estimate the distances to 
Mira-like variables through PMR after considering the reduction in 
their apparent magnitude by the circumstellar extinction. Figure [4] 
is a plot to show that there is a strong correlation between the ob- 
served near-infrared color and the deviation. This result suggests 
that the near-infrared color can be a good indicator to probe the cir- 
cumstellar extinction of Miras. To formulate the relation between 
the near-infrared colors and the deviation, we calculate the mean 
deviation (< d t >) of Miras in the x = 0.05 mag bin and its cor- 
responding standard deviation (cr^.), where the suffix i stands for 
the ;-th bin and x represents the employed color in magnitude. The 
cyan points with error bars in Figure [4] show the calculated < rf, > 
and (Tdi • Then we fit a quadratic function to the cyan points in the 
form of d = ax 2 + bx + c by using Levenberg-Marquardt algorithm 
with giving statistical weights proportional to 1/crJ.. The fitting re- 
sults are shown in the solid line in Figure [4] and the calculated 
coefficients are summarized in Table|4] 

The amount of deviation itself has important physical mean- 
ings. It can be used to roughly estimate the mass-loss rate by mak- 
ing a further assumption that the standard interstellar extinction law 
can be equally applied to Mira's circumstellar environment. This 
assumption should be inappropriate especially for C-rich Miras, be- 
cause the interstellar environment is O-rich in the first place. Also, 
the extinction law in the circumstellar environment looks indeed 
different from that in the interstellar environment, to the extent that 
the Figure [5] shows (see next). At the level of order of magnitude 
estimate, however, we suppose the assumption is valid. The ob- 
served K s band magnitude of most of the infrared Miras are more 
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Figure 5. The correlation between the pulsation period and the ratio of devi- 
ations at two different wavebands. Black and Red dots represent O-rich and 
C-rich Miras, respectively. The horiz ontal lines show the values expecte d 
by standard interstellar extinction law ICardelli. Clayton. & M athis 1989). 

than 2 magnitude fainter than the one predicted by the period-A", 
magnitude relation. There are some extreme infrared Miras that 
are more than 4 magnitude fainter than the prediction. Using the 
standard interstellar extinction law, the 2(4) magnitude extinction 
in K s band correspond to Ay of about 17(34) mag, a nd the op- 
tical d epth at visual wavelength of about 15.7(31.3). Ivan Loor] 
d2007t) related the total mass loss rate and the optical extinction 
as, M[M e yr~ l ] = 1.5 x W-\Z/Z e )-°- 5 (L/ L e ) 0JS A° w 15 . Substituting 
typical LMC values of Z = 0.4Z o and L = 8000L o , me circumstel- 
lar extinction of A v = 15.7(31.3) mag gives M ~ 1.6(2.7) x 1(T 5 
M yr~ l . If we take L = 4000L o , instead of L = 8000L o , the mass 
loss rate is reduced by a factor of (4000/8000)° 75 ~ 0.6. At any 
rate, it is likely that many of the infrared Miras in the LMC are 
losing mass of the order of about 10~ 5 M B yr~ l , comparable to the 
Galactic (i.e., presumably more metal-rich) Miras. Interestingly, 
this rather rough e s timati on of mass-loss rate is consistent with 
iGroenewegen et all d2009l) . who fit model spectral energy distri- 
bution to the observed mid-infrared spectra to calculate mass-loss 
rates of the same stars. 

The ratios of deviations in any two different wavebands also 
give us insights on the properties (e.g., grain size and composition) 
of the circumstellar dust. The basic idea is like the interstellar ex- 
tinction law that varies with the properties of interstellar dust. Fig- 
ure|5]is a plot to show the relation between the ratios of deviations 
and pulsation period of Miras. For c omparison purpose, the ratio 
of sta ndard interstellar extinction law JCardelli. Clayton, & Mathisl 
1 1989b is also shown by the horizontal line. It is clear that the ex- 
tinction law in the interstellar and circumstellar environment is dif- 
ferent. It is also clear that the ratio of deviations (hence dust prop- 
erties) changes with increasing pulsation period of C-rich Miras. 
Interestingly enough, the interstellar and circumstellar extinction 
law differs considerably for shorter period Miras, but the difference 
is little for longer period ones. It would be valuable to compare this 
observational results with the dust extinction models to quantita- 
tively constrain the circumstellar dust properties. 



4 CONCLUSIONS 

We derived the near- to mid-infrared period-magnitude and also 
period-bolometric luminosity relations of OGLE-III Mira-like vari- 



ables in the LMC. Mass-losing Mira-like variables do not fall on the 
derived relations, due to the circumstellar extinction. We showed 
that this fact can be used as an proxy to study the properties of the 
circumstellar dust, and also to estimate mass-loss rate. We found 
that the dust properties change with increasing pulsation period. 
A comparison between the observational facts and dust extinction 
model would reveal how dust grain evolves during the Mira phase. 
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